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Abstract 
(N 

0^ ■ We consider possible violations of the equivalence principle through 

the exchange of a light 'dilaton- like' scalar field. Using recent work 
. i on the quark-mass dependence of nuclear binding, we find that the 

^*" ' dilaton-quark-mass coupling induces significant equivalence-principle- 



o 
o 



violating effects varying like the inverse cubic root of the atomic num- 
ber - A -1 / 3 . We provide a general parameterization of the scalar 
couplings, but argue that two parameters are likely to dominate the 
equivalence-principle phenomenology. We indicate the implications of 
S^ . this framework for comparing the sensitivities of current and planned 

!h ' experimental tests of the equivalence principle. 

a ■ 



1 Introduction 

At the heart of the theory of General Relativity is Einstein's Equivalence 
Principle (EP). The weak Equivalence Principle predicts the composition 
independence of the accelerations of test masses in a gravitational field. 
This has been probed at a present sensitivity of 

^ ~ io- 13 (1) 

a 

in innovative and difficult experiments [HE]- Further tests of this principle 
remain important and relevant for new physics O H] . We are fortunate that 
there are several initiatives to push the sensitivity several orders of mag- 
nitude further using new space-based experiments such as MICROSCOPE 
[5], the Galileo Galilei project [6] and STEP [7j as well as new types of 
experiments using cold atoms 0(5] and sub-orbital rockets [TO] . 

One possible source of EP violation is a very liglro scalar field with 
a coupling to matter that is weaker than gravitational strength. We will 
refer to these generically as 'dilatons', although they may have origins other 
than string theory or models involving dilation symmetry. As will become 
clear below, we will phenomenologically define a 'dilaton' as a scalar field <j) 
whose couplings to matter effectively introduce a <j> dependence in the basic 
dimensionless constants of Nature (such as the fine-structure constant etc.). 
String theory may have such scalars in the low energy limit ( string dilaton, 
moduli), and these can naturally lead to EP violation at a sizeable level 
[HI E21 [TU [TU [15] . Likewise, theories of quintessence predict a light scalar, 
as do theories with continuously varying coupling constants as well as some 
theories of dark matter. While scalars lead to an attractive interaction, like 
usual gravity, they do not couple universally to all forms of energy in the 
same way as in general relativity. Thus we expect differences in the forces 
for different elements. 

Additionally, independently of any specific theoretical model one might 
argue (along the 'anthropic' approach to the issue of a possibly extremely 
vast 'multiverse' of cosmological and/or string backgrounds) that: (i) the 
'Equivalence Principle' is not a fundamental symmetry principle of Nature 
(e.g. it is 'violated' in any theory containing very light scalars); (ii) the level 
rj ~ Aa/a of EP violation can be expected to vary, quasi randomly, within 
some range of order unity, over the full multiverse of possible (cosmological 
and/or theoretical) backgrounds; (hi) as there is probably a maximal level 



We will generally assume in the following that the scalar field we consider is essentially 
massless on the scales that we discuss. 



of EP violation, say < 77* <C 1, which is compatible with the develop- 
ment of life (and of physicists worrying about the EP), one should a priori 
expect to observe, in our local environment, an EP violation 77 of order of 
77*. It is a challenge to give a precise estimate (or at least upper bound) 
of 77*. We note, however, that this is a scientifically rather well-posed chal- 
lenge. For instance, one of the necessary conditions for the existence of life 
is the existence of solar-like planetary systems stable over billions of years. 
A sufficiently large 77 7^ will jeopardize this stability, notably under the 
influence of external, passing stars. The current very small level of EP vio- 
lation ensures that stars passing at a distance D disturb the inner dynamics 
of the solar system only through tidal effects that decrease like D~ 5 . An 
EP violation 77 would increase this disturbing effect to a level oc r]R~ 2 . It is 
also a well-posed question to determine the level 77 which would destabilize 
the solar system through internal EP-violating gravitational effects. 

Independently of these various motivations, our work here will discuss 
the general type of composition-dependence of EP violation that is entailed 
by the existence of a light dilaton-like field. The theoretical challenge is to 
connect the basic couplings of the dilaton Lagrangian to the properties of 
real atomic systems. 

Our work starts in Section 2 with a review of EP violations, and a general 
parameterization of possible dilaton couplings, Eq. (fT2jh Section 3 connects 
dilaton coupling parameters with the other couplings of the Standard Model, 
which is preparation for understanding the effects of the dilaton couplings. 
Section 4 is our analysis of the effects in nuclear binding, while Section 5 is 
a summary of the effects within individual nucleons, and Section 6 describes 
electromagnetic effects. In Section 7, we collect the results of the previous 
sections and give a complete treatment of the phenomenology of equiva- 
lence principle violations, including comparisons with existing experiments. 
Section 8 provides a guide to experimental sensitivities for existing and fu- 
ture experiments. Experimenters who are willing to forgo the theoretical 
development of Section 3-6 can go directly to Section 7-8 or can consult our 
shorter paper [16] in which we have collected our most phenomenologically 
useful results. In particular, Section 7.3 contains what is probably the most 
useful parameterization of our results and Section 7.4 discusses the present 
experimental constraints. Section 9 is a brief summary. 



2 Formalism 

2.1 EP violation 

Let us start by recalling that a massless dilaton modifies the Newtonian 
interaction between a mass A and a mass B, into the form (see, e.g. |14| ) 

If the dilaton mass is important the second term includes an extra exponen- 
tial factor exp(— m^TAB)- In this interaction potential, the scalar coupling 
to matter is measured by the dimensionless factor 

1 d[Km A ((f>)} 

a A = ^i xl • ( 3 

K z m,A 0(p 

Here, k = \/4vrG is the inverse of the Planck masqj so that the product Km a 
is dimensionless. This ensures that this definition of a a is valid in any choice 
of units, even if these units are such that k depends on eft (as in the so-called 
'string frame'). In the following, we shall generally assume that we work in 
the 'Einstein frame' where the (bare) Newton constant G is independent of 
<fi. The above expression for the dimensionless scalar coupling a. a has been 
written in terms of a canonically normalized scalar field, with kinetic term 
[using the signature (+, — , — , — )] 

c 4> = \(d<t ) f + --- (4) 

Evidently, a small mass term for the dilaton can readily be added if desired. 
It can also be convenient to work with the dimensionless scalar field 

if = K(f), (5) 

whose kinetic term is related to the Einstein-Hilbert action via 

When using <p the definition of the dimensionless scalar coupling reads 

aA = 8j • (7) 



We use units such that c — 1 — H. 



In terms of the a^'s, the violation of the (weak) EP , i.e. the fractional 
difference between the accelerations of two bodies A and B falling in the 
gravitational field generated by an external body E, reads 

Aa\ n a A -a B {a A -a B )a E , . , , 

= 2 ■ = =- ~ (a a ~ a B )oLE. (8) 

a J ab a A + a B 1 + ji{aA + a B )a E 

In the last (approximate) equation we have assumed that the a's are small, 
so that one can neglect the term \(ola + ct B )a E in the denominator. 

Our aim here is to provide a general analysis of the possible EP viola- 
tions in experiments comparing the free fall accelerations of atoms (and/or 
nuclei). Most of the effort needed for such an analysis is now understood 
P2J [131 HS1 H7] , and we will use it below. However, one aspect of this analysis 
has been far less well-studied and understood, namely the contribution to 
EP violation coming from the possible ^-dependence of the nuclear binding 
energy. The aim of this paper will mainly be to assess the form of this 
contribution, coming from the quark mass contribution to nuclear binding^. 
Actually, our conclusion will be that this contribution is, possibly in competi- 
tion with Coulomb-binding effects, likely to dominate the atom-dependence 
of the EP violation signal (jSJ) • 

To motivate our general analysis, let us start by noting that the mass of 
an atom can be decomposed as 

m(Atom A ) = m A = m r A est mass + £ binding (9) 

where 

m rcst mass = ^ + ^^ + ^ ^ 

is the rest-mass contribution to the mass of an atom [Z denoting the atomic 
number and iV the number of neutrons), and where E hmdmg is the binding 
energy of the atom, which is dominated by the binding energy of the nu- 
cleus. E hmding = E% + Ei is the sum of a strong interaction contribution, 
say £3, and of an electromagnetic one, say E\ (which is dominated by the 
electromagnetic effect within the nucleus). The indices 3 and 1 are used 
here as reminders of the gauge groups underlying the considered interac- 
tions: namely, SU(3) and U(l). Note that the index A in tcia is used here 
(like in the definition of the scalar coupling a a) as a label for distinguishing 
several different atoms. It should not be confused with the mass number (or 
nucleon number) A = Z + N which we shall use below. 



3 Damour [3] and Dent [17] have highlighted this need for the study of the nuclear 
binding energies. 



2.2 The general dilaton Lagrangian 

The basic organizing principle that we shall use in our discussion is to keep 
track of the effect of all the possible <j> modifications of the terms entering 
the effective action describing physics at the scale of nuclei in their ground 
states. We have in mind here an energy scale /j, ~ 1 GeV. At such a scale, 
one has integrated out not only the effect of weak interactions, but also the 
heavy quarks c, b and t. The issue of the possible <j) sensitivity of effects 
linked to the strange quark s is more delicate. In the Appendix we argue 
that the possible EP violations linked to the <f> couplings to s are expected to 
be quite small. In the bulk of the text we shall therefore ignore s (assuming 
that its effect is taken into account by changing some of the quantities we 
discuss, notably the QCD energy scale A3). 

In this approximation, we are therefore talking about an effective action 
containing, as real particles, the electron e, the u quark, and the d quark, 
with interactions mediated by the electromagnetic (A^) and gluonic (A A ) 
fields. [Here we shall use a rescaled U(l) gauge potential, which incorpo- 
rates the electron charge e, but an unrescaled gluonic field, which does not 
incorporate the SU(3) gauge coupling #3.] Then each of the five terms in 
this effective action, say 

CcS = ~h F ^ U ~ \ F ^ FA ^ + T, [i^(A,93A A )^-miAA] , 

i=e,u,d 

(11) 

(where D(A) denotes the Dirac operator coupled to the gauge field(s) A) 
can couple to ip = K<j) with a (dimensionless) coefficient. [We assume that 
we work in the Einstein frame, with the gravity and <ft kinetic terms dis- 
played above.] This introduces five dimensionless dilaton coupling coeffi- 
cients, say d e ,d g for the couplings to the electromagnetic and gluonic field 
terms, and d me , d mu , d md for the couplings to the fermionic mass termsj. 
We shall normalize these five dimensionless dilaton coupling coefficients 
d e ,d g , d nie , d niu , d md so that they correspond (when considering the linear 
couplings to 4>) to the following interaction terms 



£int0 — K<) 



i=e,u,d 



(12) 



4 We are using here the fact that a (^—dependent coupling to the kinetic term of a 
fermion, f^ipilfiip, can be absorbed in a suitable 0— dependent rescaling of tp. 



We shall explain below the notation and our choice of normalization for 
these interaction terms. 

There are two equivalent ways of thinking about the computation of the 
scalar-matter coupling a a, Eq. ([3]). One way is to think that it is given by 
the matrix element (in the quantum state of an atom) of the operator in 
the quantum Hamiltonian (associated to the interaction Lagrangian above) 
which is linear in <p. A second way is to think that it is obtained by the 
chain rule as 

aA = dj = ^ dk a W (13) 

where KmA(k a ) is the expression of the dimensionless mass ratio Km a = 
rriA /^-Planck as a function of the dimensionless coupling constants of Na- 
ture, say k a = fci, k>2, ■ ■ ■ , &20> entering the Standard Model. Actually, be- 
cause of the limited number of terms entering the relevant low-energy action 
(jlip . there are only five relevant dimensionless constants of Nature k a cor- 
responding to the five terms in (|lip . As we shall see in detail below, the five 
terms in the interaction terms (|12p precisely correspond to introducing a <j) 
dependence in the five following dimensionless constants of Nature, 

a, KA3, nm e , Km u , Krrid, (14) 

where a = e 2 /(47r) is the fine-structure constant, A3 the QCD energy 
scale, m e the electron (pole) mass, and where m u and arid denote some 
renormalization-group-invariant measures of the light quark masses (say, 
the //-running masses taken at the multiple of A3 which is equal to 1 GeV) . 
In the next Section we shall relate our normalisation of the five dimen- 
sionless dilaton coupling parameters d a entering f|12j) to the constants (|14p . 
and explain in more detail the dependence of the mass of an atom on the 
five constants (|14p . and thereby on the five dilaton parameters d a . 

3 Relation between the dilaton coupling parame- 
ters d a and the 'constants of Nature'. 



By comparing the (^-interaction Lagrangian (]12p to the other terms in the 
effective action (|lip , we see that the meaning of the dilaton coupling coeffi- 
cients d a = d e ,dg,d me ,d mu ,d md seems clear for four of them. [Actually, we 
shall see below that the meaning of the quark-mass couplings d mi is more 
subtle, because of the renormalization group running of the quark masses, 



which is associated with the 7 mi rfg term in f)12[) .] First, the coupling d e to 
the electromagnetic field modifies the Maxwell action according to 

Cem = - l ~ A de 2 ^ F, v F^ ~ - -rj—1 -, 2 F, U F^ (15) 

4e 2 p 4(1 + d e K(f))e 2 F 

where the last equality is valid at the linear level in k</> (which is the level 
at which we define the dilaton couplings here). As we work with a rescaled 
electromagnetic field (^4 here = eA us ), the only location where the electric 
charge occurs in the Lagrangian is the one explicitly shown above. This 
allows the dilaton field to be absorbed into the following <\> dependence of 
the fine-structure constant 

a{4>) = (1 + d e K(J))a = (1 + d e tp)a. (16) 



Second, comparing fll2|) to the mass terms of the electron and the light 
quarks, we see that our normalization is such that d me ,d mu ,d md introduce 
the following (ft dependence of the e, u and d masses: 

mi(4>) = (1 + d mi K(j))mi = (1 + d mi tp)mi, (i = e,u,d). (17) 

On the other hand, the terms in f|12|) that depend on our 'dilaton-gluon' 
coupling d g call for a more subtle explanation. The choice of these coupling 
terms is such that the coefficient of d g is invariant under the renormalization- 
group (RG). As the coefficient of d mi (i.e. the mass term miipiipi) is also, 
separately, RG- invariant, our choice of normalization of the coefficients in 
(|12p gives a RG-invariant meaning to both d g and the d mi 'so 

3.1 Connection with the QCD trace anomaly 

The phenomenological consequences (for the scalar coupling to hadrons) of 
the RG-invariant nature of the couplings in (|12j) can be seen in two (equiva- 
lent) ways. One way (which was used by [18] and [15]) consists in remarking 
that the definition of the d 9 -dependent terms in (|12j) is such that they cou- 
ple (f> to the anomalous part of the trace of the gluon stress-energy tensor, 
namely 

C g<t> = -d g K<t>T" (18) 



5 We are here talking about invariance under the QCD-driven running of the QCD gauge 
coupling gs, and of the masses of fermions coupled to QCD. In view of the smallness of 
the electromagnetic coupling a ~ 1/137 <C CK3, we are neglecting the RG-running driven 
by electromagnetic effects. If one wanted to take it into account, one should add to (|12[) 
additional terms linked to the QED trace anomaly. 



where [19] 

(19) 



T-ianom 
9 



h F^F^ + j^m^ 



2^3 



Here, ^(g^) = fidgz/dfj, denotes the /3 function for the running of the 
QCD coupling g^ with the (Wilsonian) sliding energy scale //, 7m (53) = 
—fj,dlnm/dfi (with a minus sign on the r.h.s.) is the (universal) anoma- 
lous dimension giving the energy-running of the masses of the QCD-coupled 
fermions, and the subscript fi at the end indicates that the operator on the 
r.h.s. must be renormalized at the running scale [x. We recall that, classi- 
cally, the trace of the gluonic stress-energy tensor vanishes (because of the 
conformal invariance of the Yang-Mills action), but that quantum effects 
linked to the necessity of regularizing the UV infinities in the product of 
gluon field strengths at the same spacetime point x introduce the (finite) 
'conformal anomaly' (]19p [19 . Then, by using the quantum version of the 
virial theoreno, one can see [181 US] that the coupling (fTBj) means that d g 
measures the coupling of (j) to the part of the total mass-energy of the consid- 
ered hadron which is due to the (renormalized) gluonic field energy, say M g 
(where M g can be defined by subtracting from the total mass both the non- 
anomalous mass-term contributions (%2i m i' l I ) i' , Pi)> an d the electromagnetic 
one). 

3.2 Renormalization group analysis 

A second way of discussing the consequences (for the scalar coupling to 
hadrons) of our normalization of couplings in (|12p is phenomenologically 
illuminating. It consists in noting that our RG-invariant definitions are 
equivalent to very simple consequences for the <j) dependences of both the 
QCD mass scale A3, and the values of the quark masses at the scale [i = 
A3. [Note that both A3 and ttij(A3) are RG-invariantly defined quantities.] 
Let us start by defining the QCD mass scale A3 as being the mass scale 
at which the running QCD coupling g^{ii) reaches some fixed, reference 
dimensionless number of order unity, say g* = 2.5. [This numerical value, 
which corresponds to a* = g^/(47r) = 0.5, is approximately reached when 
the running scale [i ~ 1 GeV (see, e.g., the figure giving a s (/i) in the QCD 
review in [20J).] This definition of A3 can be re-expressed in terms of the 
value g c = 53 (A c ) of 53 at some high-energy 'cut-off' scale A c (which could 



We recall that this theorem says that the space integral of the spatial components 
of the total stress-energy tensor T^ t = Tjf" + T^ M + T^ tter vanishes in an equilibrium 
bound state. 



be the Planck scale, or the string scale) by integrating the j3 equation giving 
the running of <&, dliifj, = dg% / ',$3(53) , so that: 

lnA 3 (A C)5c ) = lnA c - F" -^- (20) 

Jg* P3{93) 

The expression (|20p defines A3 as a function of A c and g c . If we assume 
for simplicity that the chosen cut-off A c does not depend (in the Einstein 
frame) on (ft, the result (|20p shows that A3 will inherit a (ft dependence from 
any eventual (ft dependence of g c according to (denoting j3 c = j3z{g c )) 

<91nA 3 ^ g c dlng c 
dip ~ (3 C dp ( ) 

Similarly, the integration of the RG equation for a running fermionic mass 
rrii, namely dlnrrn = —dg^ m {g^) / ^3(53) yields the following expression for 
the value of mi at the QCD scale, In 77^ (A3): 

lnm^Ag) = lnm 4 (A c ) + F° ^4^3 (22) 

Jg, P3(93) 

Differentiating this result w.r.t. <p> then shows that the logarithmic deriva- 
tive of mj(A3) w.r.t. p> is the sum of two separate contributions, namely 
(denoting >y c = -f m {g c )) 

glnmj(A 3 ) _ ainm^Ae) g c j c dlng c 
dip dp j3 c dip 

On the other hand, by comparing^ the ^-dependent terms in (|12p to the 
basic action (jlip (both being considered at the cut-off scale A c ), we see that 
the coefficients d g and d mi have the effect of adding some ^-dependence in 
the values of g c and mj(A c ) of the form 

_^£ = - d ^L _±JiL = d + lcd 24 

dp y g c dp 

Inserting these results in the 99-derivatives of A3 and mj(A3) derived above, 
finally leads (thanks to the cancellation of the 7 c -dependent contribution in 
the derivative of the masses) to the simple results 

dlnA 3 <91nmi(A 3 ) 

-Q^ = d 9 , Q(p =d mv (25) 



7 In doing this comparison it is useful, as explained above for the Maxwell action, to 
provisionally use a 'geometric' normalization of the gluon field, i.e. to absorb g-j, in A . 



Summarizing: the physical meaning of the five dilaton-coupling coeffi- 
cients d a = d e ,d g , d me , d mu , d md is (at the linear level in cp) to introduce a 
^-dependence in the parameters entering the low-energy physics of the form 



A 3 (¥>) = 


= (l + d g if)A 3 , 


a(tp) = 


= (l + d e <p)a, 


m e ((p) = 


= (l + d me (p)m e , 


r^(A 3 )](¥>) = 


= (l + d mi (p)rni(A 3 ),i 



u,d. (26) 



3.3 Ratios of dimensional parameters 



Note that a consequence of these equations is that the dimensionless 
ratios m e /As, m u (A3)/A3, md(A3)/A3 depend on cp through the ratios 
(1 + d mi ip)/(l + dgip) ~ (1 + (d mi — d g )(p). In other words, the cp sensi- 
tivity of these dimensionless ratios is 

dln[m t (A 3 )/A 3 ] 

— = d mi - d g . (27) 

Note that this involves only the differences d mi — d g . In particular, when 
the mass couplings d mi are taken to be all equal to d g , the effect of the 
cp couplings is equivalent to introducing a cp dependence only in A3 and 
a. This fact can also be seen by means of the formulation (I18j) of the 
d g coupling. Indeed, when d mi = d g the sum of ([15]) and of the mass- 
term couplings is equivalent to having a coupling between <p and the sum 
of the anomalous, T^ nom , and of the non- anomalous, J'nonanom^ p ar ^ s Q f ^^g 
trace of the total stress-energy tensor. Therefore, modulo electromagnetic 
effects, this would imply that <fi couples to the trace of the total stress-energy 
tensor, i.e. (using the virial theorem) that cp couples to the total mass of the 
hadron. In this particular case, the only violations of the EP would come 
from electromagnetic effects. 

However, in view of the fact that the physics which determines (in the 
Standard Model) the masses of the leptons and quarks involves the sym- 
metry breaking of the electroweak sector, and, in particular, the VEV of 
the Higgs field, it does not seem a priori likely that a fundamental theory 
describing the high-energy couplings of the dilaton can ensure such a uni- 
versal feature. From this point of view, one can consider our final results 
(|26p as useful general parametrizations of the low-energy dilaton couplings, 
independently of the complicated physics that might connect these param- 
eters to an eventual high-energy description of the cp couplings to the fields 
entering the basic Lagrangian. For example, heavy quarks do not enter the 

10 



field couplings ()12p . but they enter in the relation between the QCD scale A3 
(describing the physics at scales < 1 GeV) and the high-energy boundary 
conditions, A. c ,g c . Therefore, the parametrization of d g in (|26p implicitly 
takes into account the effect of heavy quarks. [Ref. [15] showed how to ex- 
plicitly take into account the effect of heavy quarks, and it is easily checked 
that their results are in agreement with the first equation in (j26j) .] 

We can use the above results to rewrite the expression of the scalar 
couplings to matter ([3]), ([7]) in a useful form. As the Planck scale 1/k does 
not directly enter physics at the QCD scale (besides its possible impact on 
determining A3 via Eq. (|20p ). we can always write the mass of an atom as 

m A = A 3 M A (-^,-±,-^,a), (28) 

A 3 A 3 A 3 

where M A is a dimensionless quantity, which is a function of the four indi- 
cated dimensionless quantities, say (for later convenience) 

(k u ,k d ,k e ,k a ) = (—^,-—,—^,a). (29) 

A 3 A 3 A 3 

Using this notation, the scalar coupling to matter Eq.([7D can be rewritten 
(when working in the Einstein frame) as 

a A = d g + a A , (30) 

where d g = g 3 is a universal (non EP-violating) contribution to a a, and 
where the EP-violating part a a is given by 

_ d\nM A _ 1 dM A 1 ^ dM A d\nk a 
ttA= dip ~ M A dip ~ ~M~ A *-i d\nk a dip ' (ol) 

a=u,d,e,a 



The logarithmic derivatives of the k a are given by Eq. f)26|) , so that we can 
write more explicitly a A as the following sum of four contributions: 



1 dM A __ 1 
UA ~ M A dp ~ M~ A 



ST (A ^ dM A , , 9M A 

a=u,d,e 



(32) 



3.4 Redefining the quark mass parameters 

In the following, we will find it convenient to work with the symmetric and 
antisymmetric combinations of the light quark masses, namely 

m = n( m d + rn u ) , 5m = (m d - m u ) (33) 

11 



Working in terms of m and 5m, means working in terms of mass terms of 
the form 

m d dd + m u uu = m(dd + uu) -\ — 5m(dd — uu) (34) 



which couple to the dilaton as 



C, 



d m m(dd + mm) H — ^-5m(dd — uu) 



(35) 



These definitions are such that, for instance, the coupling of ip to m is 
equivalent to a Hamiltonian coupling of the form, 

H = .... + (1 + d m ip)m(uu + dd), (36) 

i.e. to introducing a tp dependence in the average light quark mass of the 
type m(ip) = (1 + d m (p)m. 

The link between these new dilaton-coupling coefficients and the previous 
ones reads 

_ dlnm d md m d + d mu m u __ dlnSm d md m d - d mu m u 

dm = —^ = ; , dsm = — q = • (37) 

clip m d + m u dip m d - m u 

In term of this notation (|32|) reads 

i r dM A dM A 

a A = —[( dlfl -d g )m— + (d Sm -d 9 )5m— 

+{d me -dg)m e — \-d e a— — J. (38) 



As displayed in Eq. (|38j) , a a is naturally decomposed into a sum of four 
contributions, which are linear in the four dilaton couplings: d ma —d g , or d e . 
Another linear decomposition can also be applied to the various terms in cya' 
namely the one corresponding to the various terms in Eq. §§§. Regrouping 
some terms in these two possible linear decompositions, we shall find con- 
venient in our calculations (before coming back to the more theoretically 
rooted decomposition ([38]) ) to decompose a a into three contributions: 

a A = a A mwo - EM + a h A ind + a d / (39) 

where a^ mwoEM denotes the contribution coming from the terms linear in 
the quark and electron masses in the rest-mass contribution (JlUp to tua 
{without the electromagnetic contributions), where a^ 1 denotes the contri- 
bution coming from the nuclear binding energy E hmd in Eq. (J9j) , i.e. 

j_a(E<>'-y/A 3 ) 

A M A dip y ' v ' 

12 



and where a A e denotes the total electromagnetic contribution, coming both 
from the EM contributions to the masses of the nucleons, and from the 
nuclear Coulomb energy term E\ , which is a part of E in Eq. ([9]) . Note 
that a A e collects the terms in a.A which are proportional to the EM dilaton 
coupling d e , i.e. which come from the (p sensitivity of the fine-structure 
constant a. This is why we have added in the definition of a A md above the 
fact that one must keep a constant when computing it. As we shall see, 
the Coulomb energy term plays a special role in that it depends both on 
nuclear-binding effects, and on EM ones. As a consequence it will give two 
separate contributions: one to a A md and one to a A . 

4 Analysis of scalar couplings to the binding en- 
ergy of nuclei 

We will first focus on the scalar coupling to the nuclear binding energy, Eq. 
(|40p . because this term has not yet received a satisfactory treatment in the 
literature. 

When dealing with nuclear binding it is convenient to work with the 
(half) sum and difference of the light quarl^] masses, m and 6m, as in- 
troduced above. Indeed, the quark-mass dependence of nuclear binding is 
dominated by its dependence on the average light quark mass m because 
pion exchanges yield the dominant contribution to nuclear binding , and 
pion masses are proportional to m, while they are insensitive to the differ- 
ence in quark masses. [ The quark mass difference 6m is important for the 
neutron and proton masses, and will enter the computation below of the 
rest-mass contribution to EP violation.] 

As explained above, the <p dependence of dfh implies the following result 
for the 'nuclear binding energy' contribution, Eq. (|40p . to EP violation: 

-bind = (^-^ gignd 

m,A om 

In QCD, because the pion is almost a Goldstone boson of the dynamically 
broken chiral symmetry, the pion mass-squared is linear in the quark mass, 
m^ ~ &0A3 m, where 60 is a pure number. This relation is accurate in the 



As explained above, heavy quarks are assumed to have been integrated out from the 
theory, thereby producing a shift in the QCD scale A3, and its associated dilaton coupling 
d g . The effect of the strange quark, which is intermediate between heavy and light, is 
discussed in the appendix. 
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physical region, so that we can translate our formula into one involving the 
pion mass, 

-bind = (d^-d g ) m2 dE*°« 

itia om^ 

Our major task then translates into knowing the dependence of nuclear 
binding on the mass of the pion. 

The semi-empirical mass formula describes the binding energy ttia — 
m rest mass through the following terms: 

m A - m r A estmass = E h[nd , (43) 

where the nuclear binding energy is approximately described as 

E »" = -a„A + aj* + aJA^l + ac ^|_I> _ ,jfc . (44) 

The various contributions to the nuclear binding energjo are called, re- 
spectively the volume energy, the surface energy the asymmetry energy, the 
Coulomb energy and the pairing energy. [In the latter, S = \ [(— ) N + (— ) ], 
i.e. 5 = +1 for even-even nuclei, 5 = — 1 for odd-odd nuclei and 5 = oth- 
erwise.] Typical fit values for these parameters are [23] a v = 16 MeV, a s = 
17 MeV,a a = 23 MeV,a p = 12 MeV,a c = 0.717 MeV. Note that, here 
and in the following, the unit of 1 MeV is supposed to represent a fixed 
fraction of the QCD mass scale, say ~ 10~ 3 A3 if we use, as indicated above, 
a reference value g* for g^ such that A3 ~ 1 GeV.] 

The m sensitivity of _E bind comes from the rh sensitivity of the various 
coefficients a v ,a s ,a a ,a c ,a p (taken in units of A3). We shall discuss succes- 
sively the m sensitivities of: (i) a v and a s , (ii) a a , and (iii) a c . Concerning 
the pairing interaction term a p we found that it was subdominant in our 
final results because it is down by a factor of A 7 '® compared to our primary 
A dependence. Even when allowing for variations with quark mass compa- 
rable to that of the asymmetry energy we found that it is negligible in the 
end, so we drop it at this stage. 

4.1 The central nuclear force terms: a v and a s 

Let us first consider the terms proportional to a v and a s . They come from the 
isospin symmetric central nuclear force, which is the dominant contribution 



9 Please be aware of a dual notation in that the letter A is used both as a label for a 
certain type of atom, and, in the semi-empirical mass formula, as a notation for the mass 
number A = Z + N. 
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in the binding of heavy nuclei. Our previous work [2H [22] shows that this 
component has an enhanced dependence on the quark masses and hence 
it has an enhanced coupling to a dilaton. This large dependence comes 
because the central potential involves competing effects of an intermediate 
range attractive force and a shorter range repulsive force. The cancelation 
between these two effects (which are individually of order ±100 MeV per 
nucleon) lead to a binding energy which is quite small on the QCD scale 
(namely of order —10 MeV per nucleon). However, the attractive force is 
far more sensitive to pion masses because it involves two pion exchange. 
Changing the pion mass a modest amount upsets the cancelation of the two 
components and leads to a larger effect than might naively be expected. 

The central force is parameterized by two terms denoting the volume 
energy and the surface energy, 

E hind = -a v A + a s A 2/3 + residual terms (45) 

The central potential is isospin symmetric, and can involve exchanges which 
carry angular momentum quantum numbers or 1. The work of Ref. [21] 
uses a general basis of contact interactions [25J to quantify these contri- 
butions to nuclear binding. This parameterization only assumes that the 
interactions have a range which is smaller than the momentum in nuclei 
k ~ 200 MeV. The dominant contact interactions are found to be those of 
an attractive scalar and a repulsive vector, describing the integrated effects 
of the potentials. They are parameterized by strengths Gs (GV) for the 
scalar (vector) channel. We can then use the results of Ref. [24j to give the 
binding energy as a function of these strengths, normalized to their physical 
values, by defining 

Vs = jr-\ > nv = yr-\ ( 46 ) 

"S I physical ^V\ physical 

This results in 

^bind = _(i2QA - 97 A 2 / 3 )r] S + (67 A - 57A 2 / 3 )rj v + residual terms (47) 

where the numbers are in units of MeV. One can see here the cancelation 
between the primary terms as each is larger than their sum. Of these two 
contributions, our calculations indicate that it is the scalar channel (jis) 
that has the most important effect. This is because the scalar channel is 
dominated by the exchange of two pions, which is highly sensitive to the 
pion mass. While the two pion contribution is often parameterized by an 
effective sigma meson, the low energy exchange of two pions is required in 
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chiral perturbation theory and is calculablqlj. This accounts for much of 
the strength typically ascribed to the sigma |26j . The vector interaction has 
a very small low energy contribution from three pions, and estimates of the 
quark mass dependence of the mass and couplings of a massive vector boson 
indicate a tiny residual contribution [2T] . 

With these results we have argued that the main contribution is the 
variation of the scalar strength with quark mass, 



a h 2 nd 



(dfh ~ d g ) 



m A 



(12(L4 - 97 A 2 / 3 ) 



m 



w dm! 



(48) 



We use the result of Ref. [21 J, displayed in Fig Q] showing the scalar strength 
as a function of the pion mass. This variation arose almost entirely from the 
threshold modification in the two pion effects at low energy, where the chiral 
techniques are most reliable and where we expect the greatest sensitivity to 
a change in the mass [28} [29] . We can use this directly to obtain 



Vs 




phys 



Figure 1: The value of the scalar strength 775 as a function of the pion mass. 



m 



dm 



m 



dys 
n dml 



-0.35 ±0.10 



(49) 



The error bar comes from uncertainties in the chiral expansion. We will 
not display the error bar in subsequent formulas, but all results in the bind- 
ing energy carry this level of uncertainty. Our final result for the central 



10 Other estimates of mass dependence [27] have not explicitly taken into account this 
low energy effect. 
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dependence in the dilaton coupling is 

^(central = ^ ~ ^ (42A - UA^) (MeV) 

m A 

( 0.036 \ , N 

« {d^-d g )F A / 0.045 --^3-J . (50) 

In the final line we have introduced the notation 

F A = ^^ (51) 

m A 

where m amu = 931 MeV is the atomic mass unit (i.e. the nucleon mass m^ = 
939 MeV minus the average binding energy per nucleon, ~ 8 MeV). The 
factor F = Am a , X mx/ m A remains quite close to one all over the periodic table 
(modulo O(10~ 3 )). Note that our result Eq. (|5Q|) for the light-quark-mass 
(rh) dependence is significantly larger (by a factor 2.2) than the estimate 
used by Dent [UJ. Indeed, Eq. (I50j) corresponds, say for the crucial surface 
energy, to a logarithmic sensitivity <91na s /<91nm = — 34MeV/a s = —2, 
while Ref. [IT] estimated d\n.a s /d\n.m ~ —0.9. 

4.2 The asymmetry energy term: a a 

Let us now discuss the (p sensitivity of the asymmetry energy ex a a which is, 
after the volume, surface and Coulomb terms, the fourth dominant contribu- 
tion to E . The asymmetry energy has two components. The first comes 
from the Pauli principle which requires that, when there is an excess of neu- 
trons over protons, the extra neutrons must be placed into higher energy 
states than the protons. The other component is due to the nuclear force in 
which the isospin dependent interactions create a stronger attraction for an 
neutron and proton compared to two neutrons or two protons. 

The asymmetry energy has been calculated by Serot and Walecka [30] in 
the same framework that we use in our work on nuclear matter [22]. This 
takes the form 

•^ - k ' j +^4 (52) 



6 JM, 2 + k 



2 12vr 2 



where 



U --a i 1 + 2 ^) '53) 



is the nucleon mass modified by interactions in nuclear matter (with Gs < 
so that M* < mjsf)- For isoscalar nuclear matter we have 7 = 4. In meson 
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exchange models G p = g 2 p /m 2 p and Gs = —g^/rn 2 are the vector meson 

and scalar coupling strengths. The kp dependence in the second term in 

a a comes from a calculation of the nuclear density in terms of the Fermi 

momentum kp. 

As mentioned above, our estimates indicate that the mass dependence of 

the vector meson coupling strength is weak. However, the Fermi momentum 

depends on the scalar strength, which has a sizeable mass variation. The 

Fermi momentum increases as the scalar strength increases. We calculate 

this through our work on nuclear matter in which we solve for the Fermi 

momentum as a function of the scalar strength (e.g. see Fig. 4 of |22j). More 

precisely, using our approximate analytical model, Eq. (17) of [22], with the 

values G s = -355.388 GeV~ 2 , and G s = +262.89 GeV" 2 (which entail the 

phenomenologically good values a v = 15.75 Me V and kp = 1.30 fm _1 ) we 

find that 

d In kp , ,. 

8h^ a525 <54) 

Using that dependence we find that both components of the asymmetry 

energy in Eq. (|52[) vary in the same direction with the scalar strength. The 

kinetic contribution (first term) varies with a logarithmic rate ~ 2.54, while 

the other one varies like kp, i.e. with a logarithmic rate 3 x 0.525 = 1.575. 

The combination of the two contributions then varies with a rate 

d In a n , 

Combining this variation with the logarithmic mass variation Eq. (|49p of 
the scalar strength Gs then yields 

2 dOa_ _ dOa_ 2 ^S_ _ _q go _ ig M y ,z.*\ 

*dml - dG 8 "dm* ~ a ~ [ ' 

Note that our framework shows that <91na a /<91nm ~ —0.82 is rather differ- 
ent from d\n.a s /d\arh = —2. This shows again the subtlety of quark-mass 
effects in nuclear physics. 

4.3 The Coulomb energy term: a c 

The Coulomb energy also has a dependence on the strong interaction cou- 
pling terms. Dimensionally this is because the electromagnetic coupling a 
is dimensionless, so that the overall energy scale associated with a c comes 
from the nuclear interactions. Physically, this dependence is also logical be- 
cause the Coulomb energy depends on how tightly the nucleons are packed 
together. We estimate this effect in this subsection. 
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An approximate analytic expression for the coefficient of the Coulomb 
contribution to the nuclear binding energy is a c ~ (3/5)ct/ro where vq ~ 
1.2 fm is the scaled nuclear radius: ta = r^A 1 ' 3 . Writing that the total 
baryonic number within the volume of the nucleus, i.e. / os47rr^/3 (with ps = 
jkp/(6ir 2 )) is equal to A, one gets the link /cf^o = (97r/8) 1,/3 . Therefore, 
vq varies inversely proportionally to kp, so that the above result shows that 
a c ex akp. This yields a logarithmic sensitivity of a c to variations of Gs 
with the same rate as kp itself, i.e. 0.525, as quoted above. Multiplying this 
rate by the rate —0.35 of Eq. fl35J), then yields 

dE 1 Z(Z-l) Z(Z-l) , s 

m—± = -0.184a c \ ' = -0.13 \ J MeV 57 

dm A 1 / 3 A 1 ' 3 

4.4 The complete scalar coupling to the binding energy 

Combining our partial results, we finally obtain for a^ the following sum 
«A nd = {drn-d 9 )F A (58) 



X 



0.036 (A-2Z) 2 A Z(Z-1) 

0.045 - — 75- - 0.020^ -^ - 1.42 x 10" 4 \..„ ; 

A 1 / 3 A 2 A 4 / 3 



In writing this result, we have, as above, factorized Fa = Am amu /mA- 



5 Scalar couplings to the rest mass of atoms 



In this section we study the first term on the r.h.s. of Eq. (|39h . i.e. the 
contribution to a a coming from the <p sensitivity of the rest masses of the 
low-energy constituents of atoms, namely protons, neutrons and electrons 
(a la PJ). 

In view of the expression Eq. (|34p for the mass terms of the light quarks, 
we can write the masses of the nucleons as [31] 

m p = m m + a -§ + C p a, 

m n = m N3 + a + -5 + C n a, (59) 

where 171^3 is the nucleon mass in the "chiral limit" of massless lightr 1 ! 
quarks, and where the electromagnetic contributions C p a, C n a will be ig- 
nored here and treated in the next section. The quantities a and 5 in 



11 In the present treatment, we absorb in rriN3 oc A3 the EP non violating effect of the 
strange quark; see the appendix. 
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Eq. (j59|) denote the matrix elements of the isoscalar (ex dd + uu) and isovec- 
tor (oc dd — uu) terms in a neutron state: 

a = {n\rh(dd + uu)\n) 

5 = {nKrrid — m u )(dd — uu)\n) (60) 

These combinations of the quark mass contributions to the individual nucle- 
ons are reasonably well known. The isoscalar contribution is related to the 
ttN sigma term and has the value a = 45 MeV |33| . The isovector difference 
can be obtained by SU(3) sum rules 

6 = md - "?" [ma - mJ = 3.1 MeV (61) 

m s — m 

The (p sensitivity of the rest mass contribution (without the EM contribu- 
tion) of an atom, 

m rmwo.EM = Am N3 + Aa +]-{N - Z) 5 + Z m e , (62) 

comes from the fact that a ex rh(ip), 5 oc 6m((p), and from the <p dependence 
of m e . Using our general results above, we therefore have 



Zm 



-rmwo.EM = {dfn _ d) _ + {d d 

rriA 2 tua 



+ (d me -d g ) -. (63) 

m A 

Inserting the numerical values of a, 5 and m e yields 

r A — 97 

-rmwo.EM ^ F A 0.0A8(d A - d g ) + 0.0017 (d Sm - d^ ^ 

hl A 



g) -T yj.yjKJ-L! V U, dm U,g) 

5.5 x 10~ 4 (d me - dg) -j-1 . (64) 



6 Electromagnetic effects 

In this section, we review the electromagnetic coupling, which is contained 
in the Lagrangian, i.e. the contribution 

(d e ) d e dm A f . 

a\ = a-r — (65) 

rriA oa 
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The main electromagnetic effects in the atomic masses come from the 
electromagnetic shifts in the nucleon masses and from the electromagnetic 
contribution to nuclear binding, E\. 



a 



(d e ) 



d e 

m A 



7 dm P ,/A V\ dmn I 9El 

Za-—^ + {A- Z)a— h a-^— 

oa oa oa 



(66) 



We follow Gasser and Leutwyler |31| in the estimate of the electromagnetic 
portions of the proton and neutron masses 



dm r 



a- 



C p = 0.63 MeV 



dm r . 



a>- 



(~"n 



-0.13 MeV 



(67) 



da da 

The electromagnetic binding is known from the semi-empirical mass formula 



a- 



dE 1 Z(Z - 1) 



da 



A 1 / 3 



with a c = 0.717 MeV. These combine to yield 



o 



(de) 



d e F A 



Z Z(Z-l) 

-1.4 + 8.2- + 7.7^^ 



x 10" 



(68) 



(69) 



As above, the factor Fa = Am Sljmu /mA can be replaced by one in lowest 
approximation. 

7 Implications for the Equivalence Principle 

7.1 General parameterization 

Summarizing our results, the dilaton coupling to an atom can be written as 



„ ^ i ^rmwo.EM , -bind , -(^e) 

cua = d g + a a + a. a + a a 



(70) 



where a A mwo - EM is given by Eq. flM}, a b A ind by Eq. dSHD, and a^' by Eq. 
(|69p . It will be convenient for the following to rewrite this result as 



ola = d g + a A (71) 

with the decomposition 

«A = [(dm ~ dg)Qrh + (dsm ~ dg)Qs m + (d me - d g )Q me + d e Q e ] A (72) 
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where Qk a can be thought of as the ' dilaton charge' coupled to the parameter 
k a . These are given by 



Qm — Fa 



0.093 



0.036 



0.020 



(A - 2Zf 



AV3 

Qsm = Fa 
Qm e = Fa 



A 2 



1.4 x 10" 



Z(Z-l) 



0.0017 



A-2Z 



A 



5.5 x 10 



-i 



Z 



and 



Qe = F A 



A_ 

Z Z(Z-l) 

-1.4 + 8.2- + 7.7^-^ 



x 10" 4 . 



(73) 
(74) 

(75) 
(76) 



Here, as above, the factor Fa denotes Fa = ^m a mu/ m ,4 (it can be replaced 
by one in lowest approximation). 

7.2 Relation to theoretical expectations 

Note that all the various contributions to the non-universal part a a of a a = 
d g + a a contain small numerical coefficients in front of the various basic 
dilaton couplings d g ,d e ,drh,ds m ,d me . It is therefore a priori probable that 
the composition-dependent part q.a is small compared to the composition- 
independentl 12 ! part a c A A ' = d g . 

We recall that the latter composition-independent part is, in principle, 
accessible in various experimental tests of relativistic gravity. For instance, 
in the notation of tests of post-Newtonian gravity, a c A A ' = d g , is related to 
the Eddington parameter 7 via (see, e.g., [14"] ) 



7-1 



4 



i + q 



-2di 



(77) 



The most precise current test of relativistic gravity [32] constrain (7 — l)/2 

1-7 



i.e. d 2 g at the level 



u g - 



< 10" 



(78) 



12 Actually, if we define the composition-independent part of cxa by some average over 
the composition of the bodies relevant for the considered gravity tests, a c ^' will have, 
besides d g , a contribution coming from oca, and notably from terms ~ 0.1(dm — dg) 
coming fron the QCD binding of nucleons, and the nuclear binding of nuclei. To simplify 
our discussion we shall assume that these terms are small. It is enough to replace some 
of our factors d g below by a c j^' = d* g ~ d g + 0.1(dm — d g ) + ■ ■ ■ to refine our estimates. 
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Planned improved solar-system tests might improve this limit to the 10~ 7 
level. As we are going to see, and as was pointed out by many authors 
before (see, e.g. |36j), such levels are much less constraining than the ones 
accessible by experimental tests of the EP. 

By contrast to the composition-independent tests whose signals are pro- 
portional to d 2 the EP-violation signals will all be (see Eq. ([8])) proportional 
to 

ole{ola - a b) — a CA -{aA - olb) (79) 

Therefore EP signals will involve the product of d g (or rather d* = a CA ') by 
one of the other dilaton couplings entering the a^'s, i.e. they will be pro- 
portional to a combination of terms involving the following four coefficients 

d* g {drh - d g ),d* g (ds m ~ d g ),d* g (d me - d g ) or d* g d e (80) 

This raises several issues of direct phenomenological interest: (i) can, in 
principle, EP experiments measure all four (a priori independent) parame- 
ters (j80j) ?; (ii) are there theoretical arguments suggesting that, among all 
the EP signals associated to these parameters, some of them might dominate 
over the others? 

Concerning the first question (which has also been addressed in |17|). 
let us note that if we approximate the factor Fa = Am amu /mA by one ( 
and Z(Z — 1) by Z 2 ) , the composition dependence of our general dilaton 
coupling above will vary, along the periodic table, according to 

oi A-2Z (A-2Z) 2 Z 2 . , 

"* = a ° + aUs +«2— ^— + a 3 ^ +a4 ^I73 ( 81 ) 

where the five coefficients oo , . . . , 04 are linear combinations (which are easily 
read off the results above) of the four dimensionless dilaton couplings dm — 
d g ,d$ m — d g , d me — d g ,d e . Here the constant offset ao is not measurable 13 ! in 
EP experiments. By contrast, EP experiments can, in principle, measure the 
coefficients of the four different composition-dependences associated with 
01,02,03,04. Barring some degeneracies, this means that, in principle, a 
well-devised set of ideal EP experiments could measure the four theoretical 
parameters ([80]) [see, e.g., [3] for discussions of the related optimization 
of the choice of materials in EP experiments, and |17j for an example of 
the determination of four theoretical parameters from four independent EP 
data] . 



At least in our approximation Am^mu/inA — 1. If one were to keep the small fractional 
(~ 10~ 3 ) variations of the ratio Arriamu/mA, one might measure part of the ao coefficient. 
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However, EP experiments will be more likely to detect signals associated 
with functions of A and Z that vary significantly over the periodic table. 
From this point of view, two signals, among the four ones in Eq. (|8ip . are 
likely to be more prominent: namely the ones associated to the parameters 
ai and 04. Indeed, both A^ 1 ' 3 and Z 2 A~ A > 3 vary significantly along the 
periodic table. By contrast, the quantities (A — 2Z)/A and ((A — 2Z)/A) 2 
vary only mildly. Indeed, the 'valley' of stable nuclei is located along a 
specific line in the A, Z plane which is rather close to the A = 2Z (i.e. 
N = Z) straight line. Actually, in absence of the Coulomb repulsion between 
protons, the Pauli principle would favour an equal number of protons and 
neutrons (cf. the discussion of the asymmetry energy above). The Coulomb 
effects modify this in favouring a relatively small excess of neutrons over 
protons. More precisely, the bottom of the valley of stable nuclei is around 

m 

1 A 

Stable ^ 2 1 + 0.015^/3 ^ 

Using this result we see that (2Z - A)/ A ~ (1 + O.OlSyl 2 / 3 ) -1 - 1, which is 
small and whose variation with A is reduced by the small coefficient 0.015. 

In conclusion, the two EP signals that are probably most easily mea- 
surable in Eq. ()8ip are the ones associated to A~ 1 ' 3 and Z 2 A~ A > 3 . In 
previous work on the phenomenological consequences of dilaton couplings 
[31 [12] it was suggested that the EP signal would be essentially proportional 
to Z 2 A~ 4 ' 3 , i.e. related to the Coulomb energy term ex d e in the results 
above. Our analysis of the quark-mass sensitivity of nuclear binding is now 
modifying this conclusion in suggesting that the (p dependence of atomic 
masses will contain, in addition to this Coulomb-related term, another term 
(related to the quark-mass dependence of nuclear binding), with a A^ 1 ' 3 
variation over the periodic table. 

An important issue is to know whether theoretical considerations can tell 
us a priori something about the relative order of magnitude of these Coulomb 
and nuclear terms. In order to discuss this we need to know something about 
the expected relative magnitude of d* g d e versus d*(d g — dm), i.e. the relative 
magnitude of d e versus d g — d m . We shall next argue that it is theoretically 
plausible either that d e ~ d g — d m , or that d e ~ (d g — d m )/A0. 

Indeed, we have seen above that our dilaton coefficients d g ,d m ,d e where 
respectively defined as being the logarithmic derivatives of A3, m, a. On the 
other hand, it is natural to consider (at least in string theory) that a dilaton 
couples with roughly equal strengths to the various terms in the Wilsonian 
action considered at some high-energy 'cut-off' scale A c , near the string 
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scale, i.e. probably near the Planck scale mp = 1/k ~ 3.44 x 10 18 GeV. If 
this is the case, the relative magnitudes of the low-energy dilaton couplings 
d g ,drh,d e is determined by the functional dependences that relate the low- 
energy quantities As,rh, a to basic couplings at the string, or Planck, scale. 
In the case, of the fine-structure constant, though it does run, according to 
the RG, between the IR (i.e. m e ) and the GUT or Planck scale, this run- 
ning is relatively small because of the smallness of the factor (2a/37r) which 
multiplies ln(mp/m e ). As a consequence, one expects that the low-energy 
EM dilaton coupling d e is similar to its more fundamental high-energy coun- 
terpart. [This is also related to the fact that we could neglect, in our action 
(fT2J) the EM analog of the ratio /3s(gs)/g3 (i.e. fiz{gz)/g\ with geometrically 
normalized gauge fields), because /3^M(e)/e 3 is essentially constant.] The 
situation is, however, quite different for the low-energy coupling d g to the 
gluon field energy. There are two equivalent ways of seeing it. One way 
(used in [15]) precisely consists in drawing the consequences of having a 
factor Ps(gs)/gs in front of (F ) 2 (to ensure RG invariance). When com- 
paring the matching of this factor at the Planck scale, versus its meaning 
at the low-scale A3 ~ 1 GeV, one sees that d g differs from its high-energy 
counterpart by a largish factor of order 



ff3(A c 

/? 3 (A C 



K = /h. q . ^?k (83) 



where the additional factor /h.q. takes into account the effect of the heavy 
quarks [15]. The second way (used in [12]) consists in differentiating the 
expression giving A3 in terms of high-energy boundary conditions. We have 
seen above that the definition of A3 coming from the integration of the RG- 
running equation for 53 yields equivalent results, with the same appearance 
of the largish factor g^,{K c ) / (3^{K C ) . It is easily checked that this second 
way also automatically includes the effect of heavy quarks, i.e. the factor 
/h.q. in K. Actually, this second way provides a quick way to estimate the 
order of magnitude of the factor K above. Indeed, the reason why A3 is 
herarchically smaller than A c is that solving the RG-running equation leads 
to a result of the type A3 ~ A c exp(— C/gl). Differentiating this expression 
w.r.t. ip immediately shows that the amplification factor between d g and 
the high-energy dilaton coupling 9 In g\ /dip can be written as 

K = ln(A c /A 3 ) (84) 

Using, for instance, A c ~ mp = 1/k ~ 3.44 x 10 18 GeV then yields K ~ 
ln(mp/TGeV) ~ 42.7, as in Ref. |12j . and consistently with the results of 
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|15j . for the MSSM case. [We note also that the presence of this logarithmic 
enhancement factor in the dilaton coupling was pointed out in Ref. |llj .] 

When considering the low-energy dilaton coupling to the average light 
quark mass m, the second way of computing it similarly suggests that it 
will contain a large enhancement factor ~ ln(A c /m) with respect to some 
high-energy counterpart that should a priori be comparable to dlng^/dip. 
Indeed, let us recall that the quark masses are of order m q ~ fH, where H 
is the Higgs's VEV, and / a dimensionless Yukawa coupling. As we do not 
know what is the mechanism which determines (from the UV) the scale of 
the electroweak breaking (i.e. which allows for a negative squared mass for 
the Higgs at low energies), we cannot compute the sensitivity of m q to (p. 
However, it is plausible, as indicated by the 'no-scale' models [37], that H 
is related to A c , via the RG-running of (scalar) masses, by an exponential 
factor similar to the one linking A3 to A c : more precisely, in these models 
one has H ~ exp(— C"//tf), where C is a constant of order unity, and where 
ht is the Yukawa coupling of the top quark. Then, the (^-derivative of In m q 
will also contain an enhancement factor of order ln(A c /A3), i.e. of the same 
order as the enhancement K above, but probably differing by a factor of 
order unity. 

Summarizing: it seems theoretically plausible that, starting from dilaton 
couplings which are of the same order, say d c = dlng^/dip, when considered 
at the high-energy scale A c , the low-energy coupling EM d e will remain 
d e ~ d c , while d g and the various d ma will be enhanced by factors of order 
K a ~ ln(A c /m a ) ~ 40. Notably, we can expect d g ~ Kd c , and dm ~ K'd c . 
This leaves us with the problem of estimating the difference d g — d m which 
enters in composition-dependent effects. It is formally of order ~ (K—K')d c . 
We do not know to what extent there could be a compensation between K 
and K' . If such a compensation exists, i.e. if K — K ' ~ 1, instead of 
~ 40, one will have d g — d m ~ d c ~ d e . On the other hand, if K and 
K' differ by a factor of order unity (or have a different sign), we will have 
d g — d m ~ 4Qd c 3> d e . Therefore, we can only write an approximate link 
of the type d e < d g — d m . For our discussion of the relative importance 
of various EP signals, it would be too restrictive to assume that Nature 
has chosen the case where d e is significantly smaller than d g — d m . We shall 
therefore continue our discussion under the general assumption d e ~ d g — d m . 

7.3 Simplified parameterization 

Our theoretical treatment of nuclear binding effects has given us some spe- 
cific predictions for the numerical coefficients of the various contributions to 



26 



the 'dilaton charges' Qk a - To better delineate what they imply for the phe- 
nomenology of EP experiments we shall henceforth make some further ap- 
proximations. First, we replace the overall factor Fa = Am amvi /mA by one. 
This is allowed because we shall see that the leading terms in the Qk a 's vary 
by factor of a few over the periodic table, while Fa differs from one only at 
the 10 -3 level. The second approximation consists in using the approximate 
equation (|82|) to estimate various Z-dependent terms in the dilaton charges. 
Namely, using this link, and taking into account the predicted numerical co- 
efficients in the dilaton charges, one finds that the terms 0.020(yl — 2Z) 2 /A 2 
(in Qm)-, and 0.0017(^4 — 2Z)/A (in Q$ m ), are numerically sub dominant. [We 
assume here that, e.g., ds m — d g ~ dm — d g etc.] In addition, we find that 
we can replace Z/A by 1/2 in Q m< , and Q e . After these simplifications, we 
can move some left-over composition-independent numerical coefficients out 
of the Q's, and into the general composition-independent contribution d g in 
a a- 

After these approximations, we end up with 



where 
and where 

and 



a A ^ d* g + [(drn - d g )Q' m + d e Q' e ] A (85) 

d* = d g + 0.093(d A - d g ) + 0.000274 



Q^ = +7.7xlO- 4Z(Z_1) 



A 4 /3 • 

We think that these approximate expressions capture all the potentially 
dominant EP violation effects. We illustrate the variation of these approx- 
imate dilaton charges over the periodic table by giving in Table 1 their 
values for a sample of elements. [Our table considers many of the same 
elements as Table 1 of [T7], but the crucial new information we provide are 
the numerical factors in the charges, as predicted from our results. We use 
the (non-integer) atomic weights as an approximate way of averaging 14 ! the 
result over the natural isotopic composition. ] 

The two main lessons we can draw from Eq. (I85P and the numbers in 
Table 1 are: (i) Contrary to what general phenomenological considerations 
(of the type of Eq. (|8"T|) ) could suggest, there are only two dominant EP 



14 Essentially we are using the approximation (f(A)} ~ f((A}), which is valid to first 
order for a smooth function f(A). 
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Table 1: Approximate EP-violating 'dilaton charges' for a sample of mate- 
rials. These charges are averaged over the (isotopic or chemical, for Si02) 
composition. 



Material 


A 


Z 


^m 


Q'e 


Li 


7 


3 


18.88 xl(T 3 


0.345 xl0~ 3 


Be 


9 


4 


17.40 xl(T 3 


0.494 xlO" 3 


Al 


27 


13 


12.27 xl(T 3 


1.48 xl0~ 3 


Si 


28.1 


14 


12.1 xl(T 3 


1.64 xl0~ 3 


Si0 2 






13.39 xl(T 3 


1.34 xl0~ 3 


Ti 


47.9 


22 


10.28 xl0~ 3 


2.04 xl0~ 3 


Fc 


56 


26 


9.83 xl0~ 3 


2.34 xl0~ 3 


Cu 


63.6 


29 


9.47 xl0~ 3 


2.46 xl0~ 3 


Cs 


133 


55 


7.67 xl0~ 3 


3.37 xl0~ 3 


Pt 


195.1 


78 


6.95 xl0~ 3 


4.09 xl0~ 3 



violation effects: one, Q' e , coming from the ip sensitivity of the fine-structure 
constant, and the other one, Q'^, coming from the tp sensitivity of the average 
light quark mass in nuclear binding; (ii) in spite of the seemingly small 
numerical coefficient entering the Q' e term, this term can be comparable to 
the Q'^ one for heavy elements, such as Platinum or beyond. Actually, one 
should remember that it is only the variations of the Q's over the periodic 
table which matters. From this point of view, note that the total variation 
of Q'^ between Li and Pt is ~ 10 -2 , while the corresponding total variation 
of Q' e is ~ 4 x 10~ 3 . Moreover, while the variation of Q'^ is localized around 
the light elements, that of Q' e keeps increasing for heavy elements. [Formally, 
Q' e oc Z 2 /A 4 / 3 ~ A 2 / 3 , while Q'^ oc A' 1 / 3 .] 

Summarizing: our theoretical framework suggests that there are two 
dominant 'directions' for the EP-violation signals associated to a long-range 
dilaton- like field, namely 



) =(olb- a c )a E = [DrhQ'rh + D e Q' e ] BC 

a J BC 



(89) 



where [Q]bc = Qb — Qc-> an d where the ' dilaton charges' are (approx- 
imately) given by Eq. (j73|) and Eq. ()76p . The coefficients D are given 
by 

An = d* g {drh - d g ) , D e = d* g d e (90) 
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where 

d* ~ a c± ^d g + 0.093(d A - d g ) (91) 

If we were assuming that the dilaton coupling d e is much smaller than 
d m ~ d g , we could go further and conclude (in view of the numerical re- 
sults indicated in Table I) that the signal Q' e is sub-dominant w.r.t. Q' Ai . In 
that case we would end up with a uni-dimensional EP signal proportional 
to [Q' m ]BC- 

7.4 Experimental bounds 

The fact that two types of EP signals are expected to dominate allow one to 
derive simultaneous constraints on the two dominant theoretical parameters 
D m ,D e by using only two independent sets of EP experiments. We can use 
to that effect the two current EP experiments which have reached the 10 -13 
level, namely the terrestrial EotWash experiment, and the celestial Lunar 
Laser Ranging one 

The EotWash collaboration has compared the relative acceleration of Be 
and Ti in the gravitational field of the Earth [lj. The Lunar Laser Ranging 
(LLR) experiments [2] measured the differential acceleration of the Earth 
and the Moon towards the Sun. We can use our framework to translate 
the results from these two experiments on constraints on the two theoretical 
parameters D m ,D e . 

The EotWash result concerns Be (A=9, Z=4) and Ti (A=47.9, Z=22), 
and reads 

" ^ QBe - a T *)a E arth = (0.3 ± 1.8) x 1(T 13 (92) 



a 



BcTi 



Working at the two-sigma level, i.e. (0.3 ± 3.6) x 10~ 13 , and neglecting the 
central value 0.3, the rewriting of this equation in terms of the theoretical 
parameters D m ,D e yields 

10~ 3 [-7.11An - 1.55.De] = ± 3 - 6 x 10 ~ 13 ( 93 ) 

The Lunar Laser Ranging measurement constrains the relative acceleration 
of the Earth and the Moon towards the Sun: 

— ) = (aEarth - a M oon)aSun = (-1.0 ± 1.4) X 10" 13 (94) 

a J Earth Moon 

In addition to the composition dependence of the matter in these objects, 
it has the remarkable ability to test the equivalence of the gravitational self 
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energy |34j . For dilaton models where the scalar also couples to matter, it 
is the matter couplings which will be most import anil 15 !, and we will not 
consider here the gravitational couplings. The Moon has a very similar 
composition as the Earth's mantle, which is mostly silicate (primarily silicon 
and oxygen). The composition differences between the Earth and the Moon 
come primarily from the Earth's core which is dominantly iron. 

We approximate the mantle composition as being Si02, and the Earth' 
core as being iron. In addition, we follow Ref. [36J in assigning to the core a 
relative mass of 32%. Working as above at the 2-sigma level, and rewriting 
this constraint in terms of our theoretical parameters^ yields 

0.32 x 1(T 3 [3.55D A + 1.0D e ] = ±2.8 x 1(T 13 (95) 

It is interesting to notice the origin of the various numerical coefficients 
in this equation, as well as in the corresponding EotWash one above. The 
r.h. sides feature the 10 -13 sensitivity level. The l.h. sides have coefficients of 
order a few times 10 -3 , which is typical for the differences of 'dilaton charges' 
listed in Table I. In addition, the LLR l.h.s. has an extra factor 0.32 due to 
the fact that only 32% of the Earth differs in composition from the Moon. 
Finally, we need to solve two linear equations for the two unknowns D^, D e 
and this introduces an inverse determinant which will further increase the 
result for the D's. At the end of the day, if one denotes eEot = ±3.6 x 10 -10 
and €llr = ±2.8 x 10~ 10 (i.e. the two, random two-sigma errors multiplied 
by 10 3 ) the solution for D^,D e reads 

D r n = -0.625 e Eo t - 3.0 e LLR , 
D e = 2.2e Eot + U.0e LLR . (96) 

If eEot and €llr were non-zero EP violation signals, this would give us 
the values of the dilaton parameters in terms of EP data. In the present 
situation, however, eEot and eiLR are only (independent) random errors. 
This expression then shows that the LLR error is dominating the error level 
in the final result. A LLR EP measurement should be about six times 



15 Indeed, gravitational self-energy couples to the combination r\g = 4(/3 — 1) — (7 — 1) 
of post-Newtonian parameters [34]. However, this combination is theoretically predicted 
|14j to be proportional to (1 — 7)/2 ~ a CI ' ~ d^ (see above). The fact that the gravita- 
tional self energy is a very small fraction of the total mass then allows one to neglect the 
corresponding effect. 

16 Strictly speaking one should take into account the fact that the EP signal involves 
slightly different values for the 'external' cie, namely the Earth versus the Sun. For sim- 
plicity, we use here the (justified) approximation where both are close to the composition- 
independent part d* of d g . 
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below the 10~ 13 level to contribute the same error level as a terrestrial EP 
measurement at the 10 -13 level. Adding the right-hand-sides of the previous 
expressions in quadrature, finally leads to the following (two-sigma) error 
levels on our theoretical parameters 

An = ±0.87 x 1(T 9 , D e = ±4.0 x 1(T 9 (97) 

7.5 Specific models 

As we discussed above, one expects that a string-theory dilaton (or mod- 
uli) will have low-energy couplings to matter of the general form d g ~ Kd c , 
d ma — d g ~ (K a — K)d c , and d e ~ d c , where d c is some common string-scale 
dimensionless dilaton coupling, where the enhancement factors K, K a are 
expected to be comparable and of order 40, and where d e does not contain 
any significant enhancement factor. Using the EotWash-LLR-derived con- 
straints given in the preceding section, we then conclude that the string-scale 
dilaton coupling d c is constrained to be d\ < 10~ 9 / (K\K — K^\) ~ 10~ 12 . 

There are two possible attitudes towards this very stringent constraint. 
One is to conclude that all the dilaton-like scalar fields of string theory 
that are massless at tree level must acquire, via loop effects, a large enough 
mass to make them invisible in current EP experiments (i.e. rrC, < 0.2 
mm). A second possibility (suggested in [12] ) consists in assuming that 
loop effects (which depend on the VEV of the dilaton) modify the usual 
tree-level dilaton dependence (oc exp(— 2p)) of the various terms entering 
the string-scale Lagrangian into more complicated functions of p, say Bi(p), 
such that these coupling functions reach an extremum at a special value, 
say 99* of ip. Indeed, under this assumption, Damour and Polyakov [12] 
have shown that the cosmological evolution of the universe drives the VEV 
of ip towards (p*, thereby ensuring that the string-scale dilaton coupling d c , 
which is proportional to d In g^/dp, is naturally very small: "Least Coupling 
Principle" (see also Refs. [38] and [39]). More precisely, [12] showed that, 
if the extremum is located at a finite field value (p* , cosmological edvolution 
would reduce an initial dilaton coupling dj, mt by a factor typicallM 17 ! of order 
Ft ~ 10 -9 . Taking this attracting factor into account then suggests that 
the present, late-cosmological-evolution dilaton coupling coefficients are of 
order 

d g ^d g - d ma ~ 40 d e ~ 4 x 10~ 8 4 nit (98) 



17 We assume here that the curvature parameter k of the dilaton-coupling function B((p) 
is of order one. See [12] for the k dependence of the total cosmological "attracting factor" 

F t (K). 
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If we insert this result into the EP violation deduced from our results 
above, say 

(a Be ~ aTi)«Earth ~ 7 x 10~ 3 dg(d g - d A ) (99) 

we get a rough "prediction" for the level of EP violation of the order 

(a Be - a Ti )a Earth ~ 10~ 17 (4 nit ) 2 , (100) 

where tf c mt is expected to be of order unity. We note that this result is 
compatible with the current experimental tests of the EP, but that several 
planned improved EP experiments [5j El [71 [8] will be able to probe this level 
of EP violation. 

In another version of this dilaton-cosmological-attractor mechanism, the 
attractor point tp* is located at infinity in field space ("runaway dilaton" 
model [H]). This corresponds to dilaton-dependent couplings of the form 

B i ((p) = C i + b i e- ip + ... (101) 

During the cosmological evolution, the dilaton runs towards (the strong- 
coupling limit) if = +oo, exponentially suppressing its coupling to matter. 
Studying the effect of this runaway mechanism during slow-roll inflation 
allowed Ref. p3] to relate the present value of the composition-independent 
dilaton coupling a CA ' ~ d g to the amplitude 5jj ~ 5 x 10 -5 of density 
fluctuations generated during inflation. This leads to 

d g x oca. ~ 3.2^4/^ (102) 

where n denotes the power of the inflaton \ i n f ne inflationary potential, 
^(x) x X a - F° r instance, in the case of the simplest inflationary potential 
V(x) = \ rn 'ix 2 1 i- e - n = 2, the above result leads to 

4~2.5xl0-(-|) 2 (103) 

In view of our present new results, Eq. (|99|) . on the level of EP violation 
associated to such a composition-independent coupling, this corresponds to 

1-10 ( b F 



(a Be - ari)«Earth ~ 2 x 10" iU I -^ J (104) 

This is in conflict with the current EP tests, except if one assumes that the 
combination of dimensionless parameters bp/(cb\) (which was assumed in 
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|13j to be of order unity) happens to be smaller than about 1/30. In such a 
model, one would expect to see EP violations just below the currently tested 
level. Alternatively, one might interpret the constraint from current EP tests 
as suggesting that the (effective) power of the inflaton in the inflationary 
potential V(x) is l ess than n = 2. For instance, if n ~ 0, Eq. (|1U2|) implies 
rfj~6x 10~ 17 , corresponding to Aa/a ~ 4 x 10~ 19 . 

Finally, a recent work [40J suggests the existence of couplings of a light 
scalar which are quite different from the usual string-motivated ones. In 
the model of Ref. [3D] the light scalar couples only to quark mass terms, 
through mixing with the Higgs. At tree level, the couplings are 

A . . 

dmi = o~ ( 105 ) 

where A is a very small mixing parameter and ran is the mass of the Higgs 

boson. However, integrating out the heavy (t,b,c) quarks (a la [181 I15j ) 

induces gluonic couplings 

2A 
d„ = ft (106) 

The constraint of this model can be then calculated to be 



A ^ 



Km 2 H 



< 4.0 x 10~ iU . (107) 



8 Experimental sensitivities 

It can be useful to use a well-motivated parameterized theoretical model 
as a guideline for comparing the significance, and relative sensitivities, 
of different experiments. For instance, the parametrized post-Newtonian 
framework |41| played a useful role in comparing the theoretical signifi- 
cance of various composition-independent tests of relativistic gravity. Here, 
we wish to capitalize on the better understanding, explained above, of the 
coupling of a generic dilaton-like field to nuclear binding energy to pro- 
pose such parametrized frameworks for comparing different composition- 
dependent tests of gravity. Our proposal is intended as an update, or a 
specification, of previous similar proposals (see, e.g. p21[T7]). Actually, our 
proposal is two-headed. 

On the one hand, if we make minimal assumptions, and essentially 
no approximations, we propose to parameterize EP violations by means 
of the matter coupling (|7ip . which involves five parameters. One of 
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them, d g (or more accurately d* = {a a)) measures the composition- 
independent part of the matter coupling, and can, in principle, be mea- 
sured by composition-independent gravity tests. The other four parame- 
ters, d m — dg,dsm — d g , d me — d g ,d e are associated with four different types 
of EP-violation signals, associated to the four different 'dilaton charges' 
Qm,Qsm,Q me ,Qe, defined in Eqs. ([73]), flTJ]), ([75]), and ([76]). 

On the other hand, we have pointed out that two 'directions' of EP 
violations are likely to dominate the measured signals. They correspond to 
the two charges Q m and Q e , i.e. to the two dilaton parameters dm — d g and 
d e . For brevity we shall denote the first one as 

d q = d m - d g (108) 

It measures the dilaton coupling to the ratio m/A^ of the average light-quark 
mass to the QCD scale. We recall that the second one, d e is associated to 
the ip sensitivity of the fine-structure constant a = e 2 /(Air). In the same ap- 
proximation that these charges dominate, we can simplify the expression of 
the the matter coupling a a and the corresponding charges, see Eq. ([85]) . and 
the equations following it. The latter, simplified two-EP-parameter frame- 
worHlj is quite predictive, and could be useful as a guideline for comparing 
and/or planning EP experiments. Let us briefly indicate some consequences 
of our proposals. 

8.1 Composition independent constraints 

The first useful result in the simplified "reference dilaton model" is the 
expected ratio between composition-independent effects and composition- 
dependent ones. As explained above the former are essentially measured by 
the Eddington parameter 19 ! 

1 - 7 ~ 2d 2 g (109) 



18 In all, this model contains three independent parameters: d g , d q and d e . If one 
could argue that the <p sensitivity of nm is much smaller than that of KA3 one could 
even consider a much more special one-parameter guideline model keeping only d g and 
setting to zero the various mass couplings d ma as well as d e . In such a model d q — —d g 
would be fixed in terms of d g . However, the no-scale supergravity models (and their string 
realizations) rather suggest that the d ma 's contain logarithmic amplification factors which 
are comparable to the one expected to be present in d g . 

19 Here d g should more accurately be replaced by some average {a a) = d* = d g + cd q , 
with a coefficient c ~ 0.1 depending of the average composition of the considered source 
bodies. 
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while the latter are given, say, by Eq. (|99p . Note that the numerical value 
7 x 10~ 3 in the latter equation comes from the Q'^ charge difference between 
Be and Ti. We can use instead the maximal difference of 10~ 2 corresponding 
to Be and Pt. This yields the approximate link 

Aa 9 d a 1 — 7 , 

10 t^t 1 no 

a d g 2 

Note that, assuming d q ~ d g , this differs by two orders of magnitude from the 
link Aa/a ~ 10 _4 (1 — 7)/2 estimated in |12| from considering as dominant 
the EM coupling d e instead of d q . This suggests that current EP tests 
correspond to post-Newtonian tests at the level (1 — 7)/2 ~ 10 -11 , i.e. six 
orders of magnitude below the current best post-Newtonian test, namely 
the Cassini limit Eq. (1 78 p . [Using the results derived above from combining 
Eotwash and LLR data, one actually gets a constraint at the level (1— 7)/2 ~ 
10~ 9 , where the loss of a factor 100 comes from the combination of effects 
explained above.] 

8.2 Test materials 

Concerning the comparison among the sensitivities of different EP exper- 
iments, we already gave above an example of the use of our framework 
(comparison between Eotwash and LLR). Let us also mention another il- 
lustrative example. Note that each EP comparison of a pair of materials, 
say (B,C), corresponds, within our simplified framework, to looking for a 
signal of the form D • Q_bcs where D is the two-dimensional vector of dila- 
ton couplings {D^,D e ) , and Q a two-dimensional vector of dilaton-charge 
differences (Q' A ,Q' e )BC = (Q'rh'Q'e)B ~ (Q'miQ'e)c- For instance, the cur- 
rent best Eotwash comparison concerned Be and Ti, i.e. (using Table 1) 
the 'charge' vector Q T iBe = (7.11,1.55) x 10~ 3 . By contrast, the MICRO- 
SCOPE experiment plans to use a pair Ti, Pt, which corresponds to the 
charge vector QptTi = (3.33,2.04) x 10~ 3 . We see that the two choices are 
nicely complementary in that the former (using lighter elements) gives more 
weight to the m component of the EP violation, while the latter (with heav- 
ier elements) gives approximately equal weights to the rh and e directions. 

8.3 Atomic interferometry 

Special mention should be given to the sensitivity of EP experiments based 
on atomic-interferometer techniques. For instance, Ref. [8] mentions the 
possibility of comparing two isotopes of Rubidium: ( 85 Rb, 87 Rb). In such 
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a case, we wish to warn the reader that one should not blindly use the 
formulas that we have derived above, especially the approximate ones for 
(Q'rhi Q'e)- Indeed, the approximations used to simplify the charges employed 
the average link (|82|) between Z and A. This approximation is acceptable 
if one compares elements that are distant along the periodic table, but is 
definitely invalid for isotopes of the same Z. Therefore, one should start 
from our original, non-approximated expressions for the charges. 

The use of our ('exact') dilaton charges, suggests that an EP test 
comparing ( 85 Rb, 87 Rb) would correspond, in the full four-dimensional 
space of (rh,5m,m e ,e) , to a charge vector equal to Qs7 Rb 85 Rb = 
(—3.3,3.4,-0.55,-9.2) x 1CP 5 . Note that the components of this vector 
are significantly smaller than those of the charge vectors probed by the 
other experiments. The dominant direction is along e. Note also that the 
5m direction now plays a role as significant as the m one, because, besides 
the binding energies, a crucial effect in comparing two isotopes is evidently 
a change in the number of neutrons. This also shows that such experiments 
are complementary to the usual ones, in that they probe new directions in 
theory space, though it comes at the cost of the overall sensitivity. 

The atomic interferometer proposal of |9j suggests the comparison of Li 
and 133 Cs atoms. In contrast to the Rubidium experiment, these elements 
are well separated in A, Z, and our simplified charges can be used. We 
find that this comparison is quite sensitive to the dilaton couplings with 
dilaton charge vector QcsLi = (11.2,3.02) x 10~ 3 . While the experimental 
comparison of dissimilar atoms may be more difficult than the use of related 
isotopes, the sensitivity to the dilaton couplings is much increased. 

Let us also make some further comments relevant for comparing two 
isotopes which are very close in mass. Our derivation assumed that the 
semi-empirical mass formula was an accurate representation of the binding 
energies. However, this mass formula is an average, which does not al- 
ways accurately capture local fluctuations, and notably fluctuations linked 
to varying A for a fixed Z. In addition, our derivation has neglected the 
pairing term — 5a p /A 1 ' 2 , as being subdominant. However, this term might 
become very important if one were to compare isotopes with mass numbers A 
differing by an odd integer. Indeed, in that case 5 = | [(— )+(—)] changes 
by one unit between the two isotopes, and therefore yields a full contribution 
a p to their mass difference, and thereby also to the dilaton sensitivity. Ac- 
tually, we would suggest to try to take advantage of this fact by using such 
odd-related isotopes which are likely to have an enhanced sensitivity to EP 
violations. [We are aware, however, that this proposal poses both theoreti- 
cal challenges (determining the <p sensitivity of a p ), and experimental ones 
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(as the two isotopes will have a different Fermi/Bose statistics, which might 
undermine the possibility of using accurate, Bose-Einstein-Condensation- 
based, techniques).] 

8.4 Other applications 

Let us also mention that our framework can be straightforwardly applied 
to comparing (weak) equivalence principle tests to atomic-clock tests of 
the dependence of coupling constants on the gravitational potential. The 
link between these two types of tests has been discussed by several au- 
thors [HI U2J H31 [T7] . Let us indicate how it is formulated in our notation. 
The spacetime dependence of the dilaton field is approximately of the form: 
<p(x,t) = tpo(t) + (p\ oc (x,t), where <po(t) is the cosmological value of (p, and 
where 

¥>ioc(s, t) = ~Yl a E-^~ - ~a cX U(x, t) (111) 

E TE 

gives the influence of the local matter distribution, in terms of the local grav- 
itational potential U (U > 0). In the second expression, we have used the 
approximation oe — a ' 1 ' = d*. Combining this result with our parametriza- 
tion k a ((p) = (1 + d a (p)k a (0) of the <p dependence of the various constants 
k a = m/A3,(5m/A3,m e /A3,a = e 2 /(47r), we see that the local gravitational 
potential influences the values of the constants k a measured, say, on the 
Earth, according to 

k l a oc = (1 - D a U)k a ( m (t)) (H2) 

where the coefficients D a = d a a CA ' = d a d* g , i.e. D m = d* (dm — d g ) ,... ,D e = 
d* d e are the same dilaton coefficients that entered our discussion above of 
the EP tests. Then, to compute the effect of the seasonally varying U 
on, say, the frequencies of atomic clocks, one needs to know the sensitivity 
of these frequencies to variations in the k a 's (see |44j). In particular, the 
D e = ±4 x 10~ 9 two-sigma bound derived above on D e , combined with the 
yearly variation AC/ ~3x 10 -10 linked with the Earth's eccentricity, shows 
that EP tests constrain the yearly variation of the fine-structure constant 
on the Earth to be smaller than 1.2 x 10~ 18 (two sigma). This is about 40 
times smaller than the current best atomic-clock experimental sensitivity to 
the variation of a |45| . Note, however, that clock-comparison experiments 
are sensitive to different combinations of the parameters D a than EP tests 
[42] . We shall not discuss here the cosmological aspects of the variation of 
constants, which are more model-dependent. For instance, in the context of 
the dilaton-runaway model, one can relate the present rate of variation of 
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the 'constants' to the (square root of the) EP violation level, see Eq. (3.25) 
of [13]. 

Let us finally remark that it would be interesting to use the recent 
progress (reported in [2TJ [22] and here) about the quark-mass dependence 
of nuclear binding to try to derive a well-justified estimate of the quark- 
mass dependence of the crucial very low-energy neutron capture resonance 
E r ~ 0.1 eV = 1CT 7 MeV of 149 Sm. Indeed, the analysis of the Oklo data 
[HH STJ HH] shows that this resonance has not changed by more than about 
0.1 eV since the Oklo natural fission reactor was in activity 2 billion years 
ago. A naive use of our results, based on our finding that the bulk binding 
energy per nucleon, a v , varies with rh as Aa v ~ — 42Alnrfi/A3 MeV, sug- 
gests that Oklo data constrain the fractional variation of rh/ 'A3 over 2 billion 
years to the level Alnm/A3 < 10~ 7 /42 ~ 2.4 x 10~ 9 . Such a limit would 
be a very significant constraint on the possible cosmological evolution of the 
dilaton. However, it is not clear whether a detailed study of the specific 
(unstable) energy level corresponding to E r will confirm this sensitivity to 
m/A 3 . 

9 Conclusions 

We have provided a parametrized framework for the study of the equivalence 
principle 20 ! in models with light, dilaton-like scalar particles. Our general 
framework contains five independent parameters, and should be applicable 
to the low-energy limit of many models. The most novel aspect of our work 
was to provide an estimate of the effects of the dilaton coupling to nuclear 
binding energy. We have found that these couplings induce, as leading effect, 
equivalence principle violations varying with the mass number as A -1 ' 3 . The 
level of these EP violations is expected to be at least comparable to (and, 
for lighter elements, somewhat larger than) that associated to the Coulomb 
energy. 

We have also provided a simplified scalar model, containing three pa- 
rameters: one composition-independent parameter, and two composition- 
dependent ones. This model is expected to describe the dominant effects of 
the most general 5-parameter framework. We suggest to use it as a guide- 
line for comparing and planning EP experiments. We used it to combine 



20 Here we have limited our considerations to the weak equivalence principle (tests of 
the universality of free fall). However, our parametrized Lagrangian can also be used to 
study the effect of dilaton couplings on other aspects of the EP: such has clock-comparison 
experiments. 
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Eotvos and Lunar Laser Ranging data so as to constrain its two theoreti- 
cal composition-dependent parameters. We found that they are constrained 
at the 10 -9 level. This plausibly implies (in our model, and using some 
naturality assumption) a corresponding limit on composition-independent 
effects at about the same level, i.e. (1 — 7)/2 < 10~ 9 , which is four orders 
of magnitude below the best present composition-independent gravitational 
tests (Cassini experiment). 

In the happy future situation of several non-zero measurements of EP 
violations, one could check the consistency of our simplified model, which 
is quite predictive. If needed the other scalar couplings could readily be 
included to make sense of subleading effects modifying the simple predictions 
of this simplified model. 
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10 Appendix: The strange quark mass 

We are not able to provide a definitive calculation of how equivalence prin- 
ciple violations depend on the strange quark couplings. This is an area 
where there is no consensus and the motifs of the day change quickly. While 
we cannot solve this issue, we will here argue that the the strange quark 
dependence could be about or within the uncertainty that we are quoting. 

When quarks are heavy, they can be integrated out with the result sim- 
ply going into a modification of the gluonic coupling, d g . The u, d quarks 
are light, are directly involved in nucleon couplings and are clearly active 
dynamically in nucleon binding. The strange quark is intermediate in mass. 
Nucleons do not explicitly contain strange quarks, so their effects are sec- 
ondary. Certainly they couple to nucleons at some level through loop effects. 
Initial theoretical calculations suggested that these couplings could be quite 
large. However, increasingly theoretical and experimental developments are 
bounding these effects to be relatively small. 

Fortunately for equivalence principle violations, the leading manifesta- 
tions of the strange quark mass would not have an effect in any case. For 
example, the much debated contribution of the strange quark to the mass of 
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neutrons and protons [49j would not lead to the violation of the equivalence 
principle. This is because the effect is an isospin singlet and contributes 
equally to the neutron and proton, so that the total effect in an atomic 
state is proportional to A. This leads to a constant contribution to a a in- 
dependent of A, and no violation of the equivalence principle. Note that the 
large effects suggested for strange contributions to nucleon masses recently 
have been bounded by lattice computations to be consistent with zero [50] . 
In nuclear binding, the leading A dependent term does not violate the equiv- 
alence principle, and it is only the surface term that is relevant. Therefore 
the key feature to be estimated is the strange quark contribution to the 
surface binding energy. 

In discussing the binding energy it is easy to be led astray. For example 
early estimates used kaon loops in chiral perturbation theory to conclude 
that there was a very large effect |51j . However, it has become clear from 
dispersive work, such as our own, that the KK intermediate state enters 
above the region of validity of chiral calculations |52j . There are analytic 
studies that show that the reliable low energy portions from such loops are 
very small [53], and lattice studies have definitively shown that the chiral 
loop effects are not strongly present at such large masses |54] , 

The lightest intermediate states involving strangeness that can couple to 
nucleons are that of a KK intermediate state and also the vector p meson 
(an ss bound state, not to be confused with our notation for the dilaton). In 
dispersive treatments, both of these start at 1 GeV. The coupling of the p to 
nucleons is highly uncertain, and depends more on the assumptions made in 
a given calculation than in a unique piece of evidence in its favor. Moreover 
it is highly constrained by recent experiments [55] that show smaller than 
expected hidden strange couplings in nucleons. If we use a estimate which 
we find to be reasonable and which is within the constraints of present 
experiments [56j, the p effects are too small to be significant. 

However, KK intermediate states can contribute to the leading scalar 
interaction and may have a non-trivial effect. We expect from most models 
of the nuclear potential that most of the scalar strength comes from below 1 
GeV. The effect of MM intermediate states must decouple as the mass of the 
meson M gets large. If we estimate generously that KK intermediate states 
contributes 10 — 15% to the scalar strength, and we take a typical form factor 
to account for the high mass threshold of the form (A 2 + &m 2 K ) _1 (where A 
is is some typical form factor scale) , we would estimate the strange quark 
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mass dependence 

-£ - <-H - < - 10 - °- 15 > aw k ~ °- 07 " ° I0 (U3) 

using A 2 = m 2 . Comparison with Eq. (J49H indicates that this is comparable 
to the error bar that we assigned to that calculation. If the KK is positive 
as expected, a contribution of this size could lead to a 20-30 % increase in 
the coefficient of the leading A -1 ' 3 term in our final results. This is clearly 
a crude estimate, but we don't expect that it is grossly misleading. 
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